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The fraction of the nucleon spin that is carried by the quarks, AE, is computed in lattice QCD with dynamical 
staggered fermions. We obtain the value AS = 0.18 ± 0.02. 



1. INTRODUCTION 

The EMC measurement of the spin- 

dependent structure function of the proton, 
gi{x,Q'^), has provoked many speculations ||^ 
about the internal spin structure of the nucleon. 
The first moment of Q^), which through the 
operator product expansion is given by the proton 
matrix element of the axial vector current weighted 
by the square of the quark charges (modulo radia- 
tive corrections), was found to be 



The quantity AE is the axial baryonic charge of 
the nucleon. 



AS = 

< p, s| ['U7p75U + ^7^75^ -I- s7At75s] |p, s > 
= <P,s\j"^\p,s >, 



(4) 
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= 0.126 ±0.010 ±0.015 



with 



Aq =< p, s\q-/^j5q\p, s>, q = u,d, s. 



(1) 



(2) 



which in a naive wave function picture can be 
interpreted as the fraction of the nucleon spin 
that is carried by the quarks. For example, in a 
SU (6)-type model of the nucleon we would obtain 
AE = 1. The vanishingly small experimental value 
of AE is referred to as the spin crisis of the nucleon. 

The nucleon matrix element of the axial baryonic 
current can be written [W 



(5) 



where is the covariant spin vector of the pro- 
ton. The average value of Q2 of the EMC data 
is « 11 GeV^. If we combine the result (1) with 
information from hyperon decays, neutron /3-decay 
and the assumption of flavor SU{3), we obtain jsj 



<P,s\fJ'\p',s' > = 

uip, s)[Gi(fc2)7^75 - G2{k^)k^ToHp', s'), 

where k is the momentum transfer and 

Gi(0) = AE. (6) 

Unlike the octet current, jjj^ is not conserved due 
to the anomaly. In the chiral limit 



AE = Au + Ad + As = 0.04 ± 0.16. 



(3) 
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*Talk presented by G. Schierholz 



where Nf is the number of flavors. As a result, 
the form factor G2(fc^) does not have a (Goldstone 
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boson) pole at A:^ = 0. Writing 
< P, s\Nf^TTF,,,F^,\p', s' > 

where m^f is the nucleon mass, one then finds 

A(0) = Gi(0) = AS. (9) 

Thus, the quark spin fraction of the nucleon spin 
is given by the matrix element of the anomalous 
divergence of the axial baryonic current. 

2. LATTICE CALCULATION 

A meaningful lattice calculation of AS requires 
(i) the inclusion of dynamical fermions and (ii) a 
proper definition of TrF^^Ffj^i^ in order to account 
for the topological origin of the anomalous diver- 
gence. None of these requirements were fulfilled in 
earlier calculations 

The calculations in this paper are done on a 
16^ • 24 lattice at /3 = 5.35 and ma = 0.01 Qwith 
four flavors of dynamical staggered fermions us- 
ing the hybrid Monte Carlo algorithm. Our data 
sample consists of 85 configurations separated by 
five trajectories. These configurations were used 
in ref. [|| to compute the hadron mass spectrum. 
The lattice parameters correspond to a renormal- 
ization group invariant quark mass of m^'^^ = 35 
MeV in the MS scheme at a scale of 1 GeV. 

For TrFf^^Ffj^^ we use Liischer's definition 0, 

X |3 f . , . , . s (fz 

Jp{x+i>,ti,u)'^ x+ii.tiu'-'P^ x+u,tiv 

^ ^x,^\i'^<7-^x,^;y\ 

~ If(x+fi,fi) i^x+(i,ti9uS'^^}^ 

y ca: f) '^^ r) ^1 1 

^ '-'x+il,^l'-'P'-'x+il,^^ '~'x+il,fi^'^'^x+fi,fii 

+ //(s.m) [^x,tJ.9,^S'^,ti ^ 

^ ^x,ti^P^x,iJ,^ ^x,ti'^'^^x,fj.^]}^ 



(10) 



where P, R and S are certain parallel transporters 
extrapolated to the interior of the plaquettes, z G 
p, and faces, z G f. This expression proceeds from 
the principle bundle which is reconstructed from 
the lattice gauge field. The resulting topological 
charge. 
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assumes integer values as in the continuum. The 
major drawback of eq. (10) is that it involves a 
three-dimensional integral over the faces of the hy- 
percubes. The more elegant algorithm of Phillips 
and Stone ^ does not lead naturally to a charge 
density. It should be noted that in the presence of 
light dynamical fermions the topological suscepti- 
bility, xt =< > /y {V- volume), is not affected 
by dislocations ||,|l^ • In the large volume limit the 
minimal action for sl\Q\ = 1 configuration becomes 

iTii 



V 

Smm oc — In(ma) , 

P 



(12) 



which grows beyond any bound for fixed bare mass 
TO (in physical units). 

According to eqs. (8) and (9), AE is given by 



AE = lim iJJ 

p-^o ps 

X < p,s|— TrF^^F^^|0,s > . 



(13) 



On a finite lattice we cannot reach p = 0. We 
therefore shall evaluate AE at the smallest (non- 
zero) momentum transfer, i.e. \p\ = 27r/16. This 
is of the order of the pion mass. Choosing s/|s| = 
± p/IpI, we thus have to compute 

c{t) = ± n 

IpI ^ (14) 
< Bp(t) P± _TrF^,Fp,(x4) 5^(0) >, 

where _B, B are the baryon creation and annihila- 
tion operators, P± is the spin projection operator 
and 



^In some cases we explicitly state the lattice spacing a in 
order to avoid confusion. 



TrF^.F^.(a:4) = ^ e'^" TrF^,F^,(a 



(15) 
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Equation (14) contains 12 independent correlation 
functions. Averaging over all of them gives f\ 

C(t) ^ AE A+ e~"^NX4~E^{t-x^) 

, . (16) 
+AEa A- (_i)te-™A^''-^A(t-^4) _^ . . . ^ 

where AEa ('tia) is the quark spin fraction (mass) 
of the A, the opposite parity partner of the nucleon 

|6), and En, a = \Jp'^ ^ "^at a ■ ^^'^ amplitudes 
are those of the ordinary correlation func- 
tion < B{t)B{Q) >. The dots in eq. (16) stand 
for contributions from higher excitations and from 
baryons backtracking around the boundary. 

A non-trivial problem is the computation of 
TrF^yF^i^ for a given gauge field configuration. 
For gauge group SU{2) we could do one inte- 
gral in TtF^^F^i, analytically which made its 
computation just feasible. In the present case of 
gauge group SU{i) we shall make use of the fact 
that the calculation can be reduced to the case 
of SU{2) by means of the so-called reduction of 
the structure group ||l^,|lj| . The argument goes as 
follows. Since TrFp,^F^i, = d^K^, we can write 
soGi(O) oc< p,s\Ko\p,s >. It turns out Q that 
^-iiTo = ^^diKoi G Z is a topological invariant. 
As in the case of the topological charge we then 
may continuously deform the SU{3) gauge field 
to a SU{2) gauge field without changing ^^Kq. 
The reason is that 7r3(5C/(3)/5'C/(2)) = 0. The 
reduction of the SU{3) link matrices to SU{2) 
link matrices is done by a coset decomposition |Q 
U{x,ijl) — uj{x, ^)U(x, where U{x,^) G SU{3), 
U{x,fi) e SU{2). In order to be able to use a 
geometric definition of TrF^^F^^ on the SU{2) 
matrices, we must make the latter as smooth 
as possible. This is done by fixing to a maxi- 
mal SU{2) gauge, which amounts to minimizing 
^{^ ~ ReTra;(a;, fj,)) ^ For the minimizing pro- 
cedure we apply a combination of Metropolis and 
overrelaxation steps. We have checked for most 
of our configurations that this results in a gauge 
invariant TtF^^F^^. 

For B we take the wall source, while for B we 
take the ordinary local baryon operator (i.e. #1 in 

■^Note that the matrix elements between A'' and A vanish in 
this case. 

^Note that this gauge is manifestly renormalizable. 
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Figure 1. The correlation function (14), averaged 
over lattice momenta p — ± 27r/16 and spins s — 
± |s|, i — 1,2,3, where Ci is the unit vector in 
z-direction, as a function of t. Also shown is a two- 
parameter fit of eq. (16) to the data, the fit interval 
being from t = A to t ~ 7. For t > 8 the data is 
too noisy to be useful. 



ref. 1^). The current is placed at X4 = 4. The re- 
sult for the correlation function (14), averaged over 
all possible momentum and spin combinations, is 
shown in fig. 1. The data is fitted by the function 
(16) with yl+, A^, 111^ and ttia taken from a fit ||^ 
of < Bf^{t)Bf^{0) >. The result is indicated by the 
solid line. It leads to 

AS = 0.18 ±0.02, (17) 

while the quark spin fraction of the A comes out 
to be AEa = 0.22 ±0.04. 

In order to check our results for systematic er- 
rors, we have repeated the calculation for x^ = 3. 
We found the same result within errors. Thus we 
can assume that a;4 = 4 is far enough away from 
the source, so that the excited states have died out. 
Since our calculation involves slowly moving (in- 
frared) modes, it is important to test the efficiency 
of the algorithm. We have found an autocorrela- 
tion time of « 25 trajectories [|ll[. This translates 
into an autocorrelation time of five configurations. 
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3. DISCUSSION 

Our result (17) lies within the errors of the ex- 
perimental value (3). One might object that our 
calculations arc not entirely realistic, as we work 
with four light quarks instead of three and the 
nucleon matrix element is computed in the chiral 
limit. But wc have reasons to believe that none 
of these approximations will have a significant ef- 
fect. As far as the flavor dependence is concerned, 
one expects AS oc y^Nj\ which leads to a small 
correction only. The quark mass and momentum 
dependence, on the other hand, is controlled by the 
7]' mass which is large compared to the pion mass. 
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